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Abstract. We develop a non-perturbative method to derive the probability distribution P(5_r) of the density 
contrast within spherical cells in the quasi-linear regime. Indeed, since this corresponds to a rare-event limit a 
steepest-descent approximation can yield asymptotically exact results. We check that this is the case for Gaussian 
initial density fluctuations, where we recover most of the results obtained by perturbative methods from a hydro- 
dynamical description. Moreover, we correct an error which was introduced in previous works for the high-density 
tail of the pdf. This feature, which appears for power-spectra with a slope n < 0, points out the limitations of 
perturbative approaches which cannot describe the pdf VISr) for (Jij ^ 3 even in the limit cr — * 0. This break-up 
does not involve shell-crossing and it is naturally explained within our framework. Thus, our approach provides a 
rigorous treatment of the quasi-linear regime, which does not rely on the hydrodynamical approximation for the 
equations of motion. Besides, it is actually simpler and more intuitive than previous methods. Our approach can 
also be applied to non-Gaussian initial conditions. 
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1. Introduction 

In standard cosmological scenarios large-scale structures 
in the universe arise from the growth of small initial den- 
sity perturbations through gravitational instability, see 
[Peebles (1980) . Besides, the amplitude of these density 
■ fluctuations usually increases at small scales, as in the 
' CDM model ( [Peebles (1982)1 ). This leads to a hierarchi- 
cal scenario of structure formation where smaller scales 
become non-linear first. Then, at large scales or at early 
times one can use a perturbative approach to describe the 
evolution of the initial fluctuations. This is usually done 
through an hydrodynamical description (e.g.. Fry (1984) 



Goroff et al. (1986)). Thus, one describes the dark mat- 



ter as a pressure-less fluid which obeys the continuity and 
Euler equations, coupled with the Poisson equation for the 
gravitational potential. However, as soon as shell-crossing 
appears this hydrodynamical description becomes inexact 
and one can no longer associate only one velocity to each 
spatial position. This implies that the perturbative series 
must diverge for hierarchical scenarios (with no small-scale 
cutoff). Nonetheless, the perturbative results obtained by 
both hydrodynamical and Boltzmann approaches are ac- 
tually identical (e.g., ^aper 1 ). 

The disadvantage of such recursive perturbative pro- 
cedures, where one computes in serial order the successive 
terms of the perturbative expansion, is that they can only 
be used for the first few order terms (e.g., up to order 



3). Indeed, the calculations become rather heavy for high- 
order terms. Hence this method cannot be used to esti- 
mate the high-order cumulants of the probability distribu- 
tion of the density field /5(x), since the cumulant of order q 
depends on the term of order (g— 1) of the perturbative ex- 
pansion. Nevertheless, for the case of Gaussian initial con- 
ditions it has been shown that one could use the structure 
of the perturbative expansion to obtain at leading order in 
the limit a ^ all cumulants of any order of the density 
contrast Sr within spherical cells (Bernardeau 1992,1994), 
where a is the rms density fluctuation. This allows one to 
get the precise shape of the probability distribution func- 
tion (pdf) V{Sii) in the quasi- linear regime. However, this 
derivation presents several shortcomings. First, it is based 
on the perturbative expansion of the density field while 
this series actually diverges for hierarchical scenarios (e.g., 
paper I). Hence the proof of the results obtained by this 
perturbative method is not complete. Second, it does not 
apply to non-Gaussian primordial density fluctuations. 

In this article, we present a non-perturbative method 
to obtain the pdf 1^(6 u) of the density contrast in the 
quasi-linear regime. It is based on a steepest-descent ap- 
proximation which yields exact results in the asymptotic 
limit (T ^ 0. Thus, it provides a rigorous justification 
of most of the previous perturbative calculations and it 
allows us to correct an error introduced in those works. 
Besides, it is actually much more intuitive. Another ad- 
vantage of our approach is that we can also study non- 
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Gaussian primordial density fluctuations, as we discuss in 
a companion paper (paper III). 

This article is organized as follows. In Sect. we recall 
the equations of motion and we introduce the generat- 
ing functions which describe the statistical properties of 
the density field. Then, in Sect.^ we describe the steepest- 
descent method which allows us to derive the pdiV{5fi) in 
the quasi-linear regime for Gaussian initial conditions. We 
also present convenient geometrical constructions of the 
relevant generating function. Finally, in Sect.H we com- 
pare our method with previous results published in the 
literature. 



2. Functional formulation 

2.1. Equations of motion 

The gravitational dynamics of a collisionless fluid is de- 
scribed by the collisionless Boltzmann equation coupled 
with the Poisson equation. Since we consider in this ar- 
ticle the quasi-linear regime it is convenient to use the 
comoving coordinate x. Then, we deflne the impulsion p 



by: 



a X 



(1) 



where a{t) is the scale factor and we note /(x, p, t) the 
distribution fmiction. Thus, /(x, p, t)d'^xd^p is the mass 
enclosed in the phase-space element d'^xd^p. Next, we de- 
fine the perturbation 6f of the distribution function by: 



/(x, p, t) = p Snip) + p (5/(x, p, t) 



where p is the mean comoving density (it is constant with 
time) and Sd is Dirac's function. Then, the density con- 
trast (5(x, t) is simply given by: 



^(x,t)-y<5/(x,p,i)dp. 



(3) 



Finally, we define the spatial Fourier transform of the field 
Sf hy: 



(5/(x,p,t) = |dke^'^-<5/(k,p,i) 
Sf{Kp,t) 



(4) 



(27r) 



dx e"*'"'' 5/(x, p, t) 



Then, as shown in paper 1 the collisionless Boltzmann 
equation can be written: 



+.l^/dk'dpW,p')^.^(k-k',p).0 



(5) 



after we used the Poisson equation to substitute for the 
gravitational potential </>. Thus, the distribution function 
5f is fully determined by cq.dq) s upplemented with initial 
conditions. As shown in paper ] , these initial conditions 
can be defined by the linear growing mode r7(k, p, t) which 



is a solution of the linearized eq.(p|). Moreover, this linear 
solution 77(k, p,t) can be derived from the linear mode of 
the hydrodynamical equations and we have: 



7j = L'+(5Lo(k)(5D(p) - i a'^D+SLoik) 



D 



k d5_ 
kP- dp 



(P) 



(6) 



where I?+ (i) is the usual linear growing mode and 5lo (k) 
is the linear density contrast today (i.e. at redshift z = 0). 
In paper I we developed a method to obtain the solu- 
tion of eq.(p|) as a perturbative expansion over the lin- 
ear mode rj. In particular, we explained that perturbative 
results obtained from the collisionless Boltzmann equa- 
tion are identical to those derived from the hydrodynami- 
cal approach. However, these perturbative series are only 



asymptotic (i.e. the series actually diverge, see paper V). 
Nevertheless, the key point is that the initial conditions 
can be defined by the linear mode ?7(k, p, t) even if the dis- 
tribution function 5f cannot be written as a perturbative 



series over 77, as shown in paper I 



2.2. Functional Z[j] 

In order to describe gravitational clustering in the universe 
we do not need to obtain the explicit solution of eq.(|5|) 
for all possible initial conditions. Indeed, since the linear 
mode 77 which sets the initial conditions is a random field 
we are only interested in the statistical properties of the 
distribution function Sf. These are fully described by the 
functional Z[j] of the test field j(x, p, t) defined by: 



(^) Z[j] = {el'^'"^^'^*^-^^) 



(7) 



where (..) expresses the average over the initial conditions. 
If we expand the exponential in eq.(R) we can also write: 



°°^ 1 f 

«— 1 y* -^ 



9=1 

X ((5/(a;i)..<5/(a;,)) 



(8) 



where we noted u) the 7-dimensional coordinate u) — 
(x, p, i). This expression clearly shows that the determi- 
nation of the functional Z[j\ is equivalent to the deriva- 
tion of multi-time correlation functions at all orders. As 
seen in eq.(||) the initial conditions ^^(x, p, t) are fully de- 
fined by the linear density contrast today (5^0 (k) since they 
can be restricted to the linear growing mode. Moreover, 
we assume in this article that these density fluctuations 
(5Lo(k) are Gaussian. This is consistent with usual scenar- 
ios of structure formation based on the simplest inflation- 
ary models, e.g. Bardecn ct al. (1983). Then, the average 
in eq.(M) is given by a Gaussian weight over the random 
held (5Lo(k). In real space x we obtain the path-integral: 



Zb-] = (DctA^i)'/' 



[d(5Lo(x)] e 



o-Sf- 



(9) 



The normalization factor ensures that Z[0] — 1, as im- 
plied by the definition (|^) (and DetA^p is the determi- 
nant of the kernel A^q). Here 5f has to be understood 
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as the distribution function which is the solution of the 
equation of motion (H) determined by the initial condi- 
tion 77(0;) defined by SloO^) ^s in eq.(^. We also used 
the short-hand notation j.Sf = J duj j{uj).Sf{u>) and 
Slo-^w-^lo = /dxidx2 (5Lo(xi).A^d(xi,X2).(5Lo(x2). 
Note that the fields j and 6f depend on the 7-dimensional 
coordinate (x, p, t) while Slo only depends on the spatial 
coordinate x. The kernel (i.e. infinite dimensional matrix) 
A^Q is the inverse of the kernel: 

Alo(xi,x2) = ((5lo(xi)<5lo(x2)) (10) 

which fully defines the statistics of the random field (5lo(x) 
since the latter is Gaussian. We define the Fourier trans- 
form of the kernel A^q by the property: 



contrast, the path- integral IM) involves a simple Gaussian 
integration over the initial conditions but the factor j.Sf 
is not explicitly known. Nevertheless, the expression (0) 
will prove to be more convenient because we shall not 
need the explicit mapping (5/[(5lo]- Indeed, as shown in 
the next sections we shall only need particular spherical 
solutions. 

2.3. Generating functions 



■fpiy) 



-ySR\ 



/l.Aio./2 = /"dkidk2 /i(ki)*.Aio(kl,k2)./2(k2) (11) 

for any real fields /i and /2. Using eq.(0) this implies: Sr — / —— S{x.). 

Jv ^ 

ALo(ki,k2) -ydxidx2 e^i^--^--^^-^^) Alo(xi,X2) (12) 

which gives: 

ALo(ki,k2) = {2nf Po(fci) <5z3(ki -k2) (13) 

where we defined the power-spectrum Po{k) of the linear 
density contrast today by: 

((5Lo(ki)(5Lo(k2)) = Po(fci) Soiki + k^). (14) 

Then, since the inverse A^^ is defined by the property: 

/ dx A^q(xi,x) Alo(x,X2) = (5d(xi - X2) (15) 

which reads in Fourier space: 

/dk A^i(ki,k) Aio(k,k2) = {27Tf Soiki - ks) (16) 



In Sect. 2.2 we introduced the functional Z[j] which pro- 
vides all statistical properties of the stochastic field 
(5/(x, p, i). However, in practice one does not need all of 
these properties oiSf. In particular, one is often mainly in- 
terested in the pdiViSu) of the density contrast 6r within 
a spherical cell of comoving radius R, volume V: 



(19) 



Without any loss of generality we can take this cell V to 
be centered on the origin x = 0. Next, it is convenient to 
introduce the generating function ?A(y) defined by: 



ddR e-y''^ V{5r). 



(20) 



Here (..) is the average over the initial conditions. The 
last equality in eq.(|2C|) defines the pdf V{Sr). Note that 
the generating function iI}l{v) associated with the pdf 
'Pl(^l._r,) of the linear density contrast is related to Vl 
as in eq.(Efl) where the integration over 5l,r now runs 
from —00 to -fcx). The pdf can be recovered from ^j{y) 
through the inverse Laplace transform: 



J-ioc 27ri 



(21) 



we can see that the kernel A^o obtained in eq. (^_3|) can be 
inverted as: 



AZo'(ki,k2) = 



1 



P^iki) 



(5n(ki -k2) 



(17) 



Moreover, we can see from eq.(0) and eq.(p^) that the 
inverse kernel A^q is positive definite since we have: 



In fact, as shown by eq.(M) one can directly obtain the 
moments of the pdf from the generating function il}{y): 

V'(y) = E^('5^)- (22) 

Gomparing eq. (pO|) with eq.(R) we see that we have ipiv) — 
Z[Jr] with: 



5l.^ 



W^L 



dk 



Po{k) 



(18) 



, , 9(x < R) ^ . , 

Jfl(x,p,t) = -y Ooit-tR) 



(23) 



where we used ^^(^k) — S^ik)* for real fields 5l(x), see 
eq.(H). Note that Alo(xi,X2) is simply the two-point cor- 
relation function ^lo(|x2 — xi|) of the linear density field 
today. Hence Aio(xi,X2) is symmetric, homogeneous and 
isotropic, as impHed by eq.(p^). 

Note that in paper I we obtained an alternative path- 



where 6(x < R) is a, top-hat with obvious notations and Ir 
is the time which we consider. Then, instead of integrating 
over the linear density contrast today 5lo as in eq.(H) we 
can equivalently integrate over the density contrast at the 
time of interest 5l (since 5l oc Slo)- This yields: 



integral expression for the functional Z[j]. It involved an 
integration over the actual distribution Sf{u)) and the as- 
sociated weight was not Gaussian: the argument of the 
exponential contained terms of order two to four over the 
field 5f. Moreover, all the terms were explicitly known. By 



i^iy) 



(DetA 



-IN 1/2 



[d5L(x)] e 



-vSrISl 



Ul-AT'-Sl 



(24) 



where we introduced the symmetric kernel: 

Al(xi,X2) = (5l(xi)5l(x2)) = a(|x2 - xil). (25) 
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Thus, Ar is the inverse operator of the two-point correla- 



tion function ^^ of the linear density field. In eq.fcj) the 
quantity (5_r[(5l] is the exact non- linear density contrast 
Sfi over the cell V which arises from the gravitational dy- 
namics of the linear density field (5l(x), as described by 
the equation of motion (H). Note that eq.(ni]) to eq.(n8h 
also apply to A^, A^^ and P{k). 



by the minimum of the action 5* in the limit a —^ for 
a fixed y. Indeed, the contributions from other states Sl 
are exponentially damped relative to this point. Moreover, 
the steepest-descent approximation becomes exact in this 
limit. Of course, this is the reason why we performed the 



change of variable y 
eqM 



y/u" in the change ip 



Lp m 



3. Steepest-descent method 

3.1. Action S[6l] 

The calculation of path-integrals such as (|23) is in general 
a rather difficult task. However, when a parameter be- 
comes very small one may try a steepest-descent approx- 
imation. Indeed, it may happen that in such a limit the 
integral in eq. ( 24 ) becomes increasingly dominated by the 
point where the argument of the exponential is maximum 
(i.e. the minimum of the "action"). See for instance any 
textbook on Quantum Field Theory for a discussion of the 
steepest-descent approximation. In this article we consider 
the quasi-linear regime. Then, the parameter which tends 
to zero is the amplitude of the linear two-point correlation 
A/,, that is the amplitude of the linear power-spectrum 
P{k) at the time of interest. 

In order to factorize the amplitude of the two-point 
correlation A^ it is convenient to define a new generating 
function (p{y, cr) by: 



i^iy) 



^ g-vlj/o" ,'^)/'^ 



(26) 



where we noted as usual a{R) the rms linear density fiuc- 
tuation in a cell of radius R: 



a^R) ^ (SU) 



f dxi dx2 

= / Ar(Xl,X2) 

Jy V V ^ ' ' 



(27) 



Here Sl.r is the linear density contrast within the spher- 
ical cell of radius R. Hereafter, we consider the shape of 
the power-spectrum to be fixed and the second argument 
of (p{y, cr) describes the dependence of the pdf on the am- 
plitude of the power-spectrum. The generating function 
ip{y,a) provides the pdf V{5r) through eq.(21). It also 
yields the cumulants (^^)c of the density contrast through: 



v{y,^) 



E 

q=2 



(-1) 



9-1 



{S'a 



Ric 



g! 



r2(9-l) 



y^ 



where we used the fact that {Sr) 
obtain: 



0. Using eq.(||) we 



-¥'(J',<t)/o- 



(DetA 



-ni/2 

L ) 



[d(5L(x)] e 



-S/a^-(R) 



where we introduced the action S 
a^R) 



S[dL] = ySR[SL] + ^ 



Sl.Aj\5l 



(29) 



(30) 



The action S[6l] is independent of the normalization of 
the linear power-spectrum P{k) since A^ oc a^. Then, 
it is clear that the path-integral in eq.(E9) is dominated 



3.2. Spherical saddle-point 

Thus, for any y we look for the point Sl where the action 
S is minimum. The condition which expresses that 5l is 
an extremum (or a saddle-point) is: 



SS 



S{SLi^)) 



— for all x. 



(31) 



where 5/^(5l(x)) is the functional derivative with respect 
to Sl at the point x. This constraint also writes: 



S{Sr) 



cj^{R) /dx' A^1(x,x')<5l(x') = 



(32) 



since A^ is symmetric. Hereafter, we consider y to be 
real. Multiplying both sides of eq.(32) by the operator 
Al/(t^(J?) we get after a change of notations (see eq.(|5|) 
for the property which characterizes inverse operators): 



<5l(x') = 



a^{R) 



dx" Ai(x',x" 



S{Sr) 
5(<5i(x"))- 



(33) 



Since Ai(xi,X2) = ^l(|x2 — Xi|) is homogeneous and 
isotropic we can look for a solution of eq.(33) which 
is spherically symmetric. For such a spherical solution, 
(5l(x') is fully defined by the overall density contrast Sl,r' 
within spherical cells of arbitrary radius R' and volume V' 
centered on the origin. Then, taking the mean of eq.(p3|) 
over a cell V we get: 



3L,_R/ 



-y 



dx' 



dx' 



„ Ai(x',x") S{Sr) 



cT^R) <5(<5l(x")) 



(34) 



Going back to the definition of derivatives we can write 
eq.dH) as: 



(28) S 



L.R' 



lim 

£^0 e 



5r 



dx' Al(x',x") 



/• dx' Ai( 

5r[5l{^")] 



(i?) 



(35) 



We can see that both arguments of the functionals 5r 
which appear in eq.(|35|) are spherically symmetric. Hence 
we only need consider the restriction of 5r[8l{^")\ to 
spherically symmetric states, which we note 8r[5l,r"]- 
Here x" and R" are dummy variables. Thus, we write 
eq.(pa) as: 



5l,w = lim — {&R Ul,r" + eAi(i?', R")] - Sr [Sl.r"]} 

e^O ell J J 



(36) 
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where we introduced the normalized kernel: 



dxi 



aHR) 

dx2 Al(xi,X2) 

V2 



(37) 



which is independent of the normalization of the power- 
spectrum and it satisfies: Al{R,R) — 1. Thus, the linear 
density profile 5l.r" must satisfy the constraint ( |36| ) for 
any radius R' . This fully defines the saddle-point 5l.r"- 
The matter enclosed within the radius R for a spherical 
linear state Sl.r" comes from a Lagrangian comoving ra- 
dius Rl- Moreover, if there is no shell-crossing the actual 
density contrast 5r only depends on the linear density 
contrast Sl,Rl because of the spherical symmetry of the 
system. Thus, for such spherically symmetric states we 
can write; 



6r = T {Sl.r.i_ 



(38) 



i?i = (1 + bR)R' 

where the function T\bh,Rj\ is given by the usual spherical 
collapse solution of the equations of motion. For instance, 
in the case of a critical density universe T[&l) is defined 
by the implicit systems (Peebles (1980)): 

9 ^« -^"^^2 



1+^('5l) 



sml 



^r >0 



and 



^1. <0 



1 + T(8l) 



2 (1 -cos( 



- [6(0- sin 0)]^/^ 



9 (sinhr; — 77)^ 



(39) 



2 (cosh 77 -1)3 



(^L =-^ [6(smh?7-77)]^/^ 



The second equation in (pq) merely expresses the con- 
servation of mass. In order to solve eq.(|3q) we need the 
change of ^r to first order in e under the transformation: 



U.fl" 



5L^fl"+eAL(i?',i?"). 



(41) 



Under this infinitesimal transformation both quantities 8r 
and i?i are modified and we write to first order in e: 



^R -^ Sii + e (5ll^ and Rl ^ Rl + ( R 



(1) 



(42) 



Substituting into the system (38) we obtain for the first- 
order term: 



S'^R^ = T' \hLM, 



R 



(1) 



d,5 



L.R" 



AR" 



Al{R',Rl) 



The only term in the r.h.s. of eq.(M) which depends on 
R' is the factor Al{R' , Rl) hence we have: 



JL,R; 



, ^5 



L,Rl 



Al{R',Rl) 



Al{Rl,Rl) 
As a consequence, we can write: 

dSL.R" ^ 1 



^Sl. 



'%iR",RL). 



dR" "'^" Ai(i?i,i?i)ai?" 

Here we note that we also have the identity: 

a^{R")^a^{R)AL{R",R"). 

Then, since Al{Ri, R2) is symmetric over Ri 
can easily see from eq.([l7|) that: 



(45) 

(46) 

(47) 
<-> R2 we 



1 



dAL 



Al{Rl,Rl) dR 



^{R'\Rl) 



1 



%iRL)- (48) 



Hence, using eq.(M) we can write eq.( 
R' = Rl as: 



a{RL) dR 

at the point 



^L.Rl 



T'[5l,rA<jHRl)/<j\R) 



l-T'[5L,Rr.]^5L,R 



L cr{RL) dR 



mi^L) 



(49) 



This fully defines the spherically symmetric saddle-point 
(5l(x). The implicit equation ( [49[ ) determines Sl.Rl while 
the radial profile of this initial state is given by eq.(Eq). 
The radius Rl is related to R by the second equation of 
(Bq). Note that this saddle-point is an exact solution of 
the coUisionless Boltzmann equation: we do not use an 
hydrodynamical description of the system. 



(40) 3,3, Density profile of the saddle-point 

Before we derive the generating function ip{y,(7) implied 
by the spherical saddle-point obtained in the previous sec- 
tion, we need to examine its radial density profile. For in- 
stance, if the density contrast were to become larger than 
unity at radii greater than Rl this would invalidate the 
previous results since these outer shells would have col- 
lapsed and relations (Bq) would no longer hold. The ra- 
dial profile is given by eq.(p5|). Let us define the Fourier 
transform F{kR) of the top-hat of radius R by: 



Rl 



r Hx 
F{kR) = / — e*"^-" = c 
Jv ^ 

Then, we obtain from eq, 
1 



sin(/i:i?) - [kR) cos{kR) 
{kRf ■ 

dl) and eq.(|l|): 



(50) 



(43) AL(i?i,i?2) 



^RlRl 



(1) _ /;(1) 
(1) 



^r'R"" 



which yields 5j^ . From eq. ( 
and we obtain: 



Slm' = -y 



we see that Sl.r' 
P[5l.,Rl]^l{R',Rl) 



-yS_ 



(1) 

R 



aHR) 
which yields for the density profile: 

j cWi FlkR') F{kRL) P(k) 



dk F{kRi) F{kR2) P{k) (51) 



TL.R' 



(>L.Rl 



J dk F{kRL)^ P{k) 



(52) 



1 - ^' [Sl,Rl] 



fl3 

sir 



dS 



L,R" 



dR" 



(44) 



To get an idea of the radial profile implied by eq.(p2D 
it is convenient to consider the case of a power-law linear 



p. Valageas: Dynamics of gravitational clustering II. Steepest-descent method for the quasi-linear regime. 



ir^ 



O ^ 



-3 



n=-l 



-1 



n=0.5 



-0.5 0.5 

Log(R'/RL) 



Fig. 1. Cumulative linear density profile Sl.r'/Sl,b.l oi 
the spherical saddle-point. The solid line corresponds to 
n — —1 and the dashed-line to n = 0.5. 



so that eq.(|^ and eq.(|^ hold for all radii R'. For illus- 
trative purposes we show in Fig.n^ the cumulative linear 
density profile of the spherical saddle-point for the cases 
n = —1 and n = 0.5. Note that for n > the amplitude of 
the density perturbation actually shows a slow decline for 
R' < Rl. Moreover, for n > the local density contrast 
(5i(x) changes sign at |x| = Rl (so that |(5i.jj/| decreases 
faster than R'~^ at large radii) while for n = the profile 
5l (x) is simply a top-hat of radius Rl ■ 

3.4. Generating function ip{y) 

The spherical saddle-point we obtained in the previous 
section provides the asymptotic behaviour of the pdf 
V{Sfj) in the limit cr — > 0. Thus, it yields the limiting 
generating function (p{y) defined by: 



power-spectrum P{k) ex /c". Then, we can write eq.(p2[) (^(y) = viv^'^ = 0) 

i?' \ -3/2 joo^^ k^-^h,2{kR')h,2{kRL) 



(60) 



as: 



TL,Rl 



Rl 



/(f dfcfc«-l J3/2(fci?i) 



■(53) 



To write ^p{y) we first need the value Sy of the action S[5l\ 
at the saddle-point, which is given by eq.(|30|). Since the 



where we used the fact that, using eq.(|50|), F{x) can also 
be written: 



saddle-point obtained in Sect. 3.2 is spherically symmetric 
we can check from cq.(E5h and the definition (p^) that it 
obeys: 



F(x) = 3 J|a;-3/2j3/2(a:). 



(54) 



(5l(x) = 5l,Rl 



1 



dx' Ai(x,x') 



Here J3/2{x) is the standard Bessel function of order 3/2. 
The expression (53) also reads: 



}L,R' 






2Rl Y+' 2Fi(^,2;4; 



ALiRL,RL)Jv^ Vl 
Substituting this expression into eq. 
1 a^iR) 



a^{R) ■ 
ly) we obtain: 



{Rl+B'V 



Sy ^vTISlml] 



2 a^RL) 



SIr. 



(61) 



(62) 



Sl,r. \Rl+R'J 2^1 (^,2; 4; 1) 



(55) 



where 2F1 is Gauss' Hypergeometric function, see 



Gradshteyn fc Ryzhik (1965)| (§6.576.2). Th en, using the 

(56) 



relations (see pradshteyn fc Ryzhik (1965)| , §9.122.1 

+ 3 



,2;4;1 = 



24 



(l-n)(3-n) 



2F1 
and 
2F1 
we obtain the asymptotic expressions: 



Then, applying the steepest-descent method we approxi- 
mate the path- integral in eq. (|2^) by the Gaussian integra- 
tion around the spherical saddle-point. This yields: 

^^^iv.a)/a- = (DetA^i)'/' (DetM^)-^/' e-^»/"'(^) (63) 

where the minimum Sy of the action S is given by eq. ( |62|) 
while the matrix My is the Hessian of the exponent: 



,2;4;0 = 1, 



i?'^0 
and 

R' ^oc 



OL,B.' 
^L,B.L 



Sl.r 



^L.Rl 



.n (l-n)(3-n) 



^„ (l-7i)(3-n) fRL 
R' 



n+3 



(57) 



(58) 



(59) 



Mj^(xi,X2) = 



52 (5/^2) 



6{dL{^,))S{SLi^2)) 



a^{R) 5(<5l(xi))<5((5l(x2)) 



AZ'(X1,X2). 



(64) 



Note that eq.(|63[) becomes exact in the limit cr ^ if 
the spherical saddle-point found in the previous sections 
is indeed the global minimum of the action. Taking the 
logarithm of eq.(|6^) we see that the contributions to f{y) 
from the determinants vanish in the limit a — > so that 
we simply get: 



Note that for all cases of cosmological interest we have 
— 3 < rt < 1. Thus, we see that at small radii 5l,r.' re- 
mains finite, of the order of Sl.Rl 1 while at large radii it 
decreases as i?'~("+3) c>c a^{R'). The refo re, the derivation 
of the spherical saddle-point in Sect. 3.2 is valid since for 
|<^L,i?i I ^ 1 the density contrast remains small at all radii 



v(2/) = Sy = y T [5l,r.l 



riSL,R.r 



where we introduced the function t{Sl) defined by: 
- Sl a{R) 



t{5l 



a[{l + T[5L]Y/^R] 



(65) 



(66) 
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We introduced a minus sign in eq. (pq) in order to recover 
the results of the perturbative hydrodynamical approach 
(see eq.(|69|) and eq.([70|) below). Next, taking the deriva- 
tive t'((5l) we can recognize the structure of the denomi- 
nator in eq.(E9h and we can write eq.(E9) as: 



(67) 



ddL 



Thus, the intersection of the curves rdr/dSL and 
—yJ^'(5L) determines the saddle-point Sl^h^^, or r, associ- 
ated with a given value of y. We can simplify this implicit 
system by introducing the function G{t): 



g{T) = T[SL{T)]^Sf 



(68) 



where Sl{t) is defined by eq. (puj) . Using eq.(|66|) we can 
see that G{t) is also defined by the implicit relation: 



G{t)^T 



a{R) 



(69) 



Then, eq.(|67|) and eq.(|65|) can be written: 



(T = -yg'{T) 



(70) 



v{y) ^yG{T) + — 



This implicit system fully defines the generating func- 
tion f{y), with Q{t) defined in eq.(p39|). As discussed in 
Sect.0, this agrees with the results obtained from the 
usual perturbative hydrodynamical approach (Bernardeau 
(1992, 1994)). Note that for a power-law power-spectrum 
P{k) DC fc" we have a{R) oc i?-(»+3)/2 go that eq.(|9|) 
simplifies to: 



q{t)=t -T{i + g[T]) 



-(n+3)/6 



(71) 



Finally, we note that the definition ( |39D for the function 
J-"((5l) breaks down for (5l ^ 1 where shell-crossing occurs. 
In fact, beyond this point the density contrast 5r also 
depends on the density profile of the initial condition so 
that we can no longer write a relation of the form (p8|). 
Hence, our results only hold for (5l ^ 1 (i.e. as long as 
there is no shell-crossing beyond the radius R). However, 
in the quasi-linear regime cr <C 1 the typical values of 5^ 
are small: \5l\ ^ cr{R). As seen from eq.(|66|) and ( [70| ) this 
corresponds to \y\ ~ \t\ ~ o-{R) <C 1. Thus, in the quasi- 
linear regime the pdf ViSp) should be well-described by 
the system ([70[). Indeed, this limit corresponds to cr ^ 
for a fixed y, which we can take to be small but finite (e.g., 
I2/UO.I). 

However, before we can reach this conclusion we must 
check two points. First, we must make sure that the 



saddle-point 5l we obtained in Sect. 3.2 is indeed a min- 



imum of the action (and not a maximum for instance). 
Second, we must ensure that it is the global minimum 
(and not a mere local minimum). 

To check the first point we simply need to make sure 
that the Hessian VF(xi,X2) of the action S[5l\ is positive 



definite at this point. Since W = a^{R)My where My is 
the kernel defined in cq.(pj) we have: 



T^(X1,X2) 



y 



s^s 



(5(<5i(xi))<5(<5L(x2)) 



S{Sl{^i))S{Sl{^2)) 



(72(i?) A-i(xi,x2). 



(72) 



(73) 



The kernel (j'^{R) A£^(xi,X2) is of order unity since it is 
independent of the normalization of the power-spectrum 
and of y. Moreover, it is positive definite, as shown by 
eq.(18). On the other hand, the first term in the r.h.s. of 
eq.(73) vanishes as y for y — > (indeed, 5l,Rl ^ for 
2/ ^ so that the second derivative tends to its value 
at the point (5l = 0). Then, if we define the determinant 
D{y) = det(W^) we have D{Q) > 0. Since D{y) is a con- 
tinuous function of y we can conclude that D{y) > over 
a finite range of y around 0, which implies that W re- 
mains positive definite over a finite range of y. Hence, for 
small y the saddle-point we obtained in Sect. 3.2 is indeed 
a minimum of the action S. 

Next, we must show that this local minimum is in fact 
a global minimum of the action. As explained above, we 
take y to be small (but finite), \y\ <C 1, since we study the 
quasi-linear regime. On the other hand, the second term in 
eq. (p3) is of order unity. Then, we see that if there exists 
another local minimum 5^ of S it must be at least of 
order unity. Indeed, in the neighbourhood of the spherical 
saddle-point 5^ where the Hessian W is dominated by 
a^A^^ there can be no other saddle-point. 

Let us first consider the case of positive y. As seen 
from eq.(|70|) and eq.(p6|) this corresponds to positive r 
and negative 5l (since Q{t) is a decreasing function of 
r). In fact, this could be directly seen from eq.(|30|) which 
clearly shows that in order to minimize the action S with 
a positive y we must have 5h < Q. Moreover, since the 
density p must be positive we have the constraint 5r > — 1 . 
Hence we obtain from eq. (pOJ) : 



\R) 



5l.A-, 



(74) 



2/>0: S[SL]>-y+- 

Since \y\ <^ 1 and 5^ 

S[S'^] ~ ia^iR)/2) {6'^.Al\S'^) <; 1 since A£" is positive 

definite. On the other hand, we can check from eq.(p5|) 



1 we see from eq.(|74|) that 
1 — 



-y /2 < for the spherical saddle-point 



that S[Sl] 

we obtained in Sectj3.2| (indeed, for \y\ ^ 1 we have: 
T ~ y ~ —5l,Rl — ~^r)- Hence if there exists another 
minimum 5^ of S it obeys S[S'^] > S[6l]- Thus, we con- 
clude that for small positive y the saddle-point obtained 
in Sect. 3.2 is the global minimum of the action, which 
justifies the previous calculation. 

Finally, we consider the case y < 0, which corresponds 
to 6b, > 0. This case is more difficult since there is no 
upper bound for Sr and for large Sl we no longer have a 
relation of the form (p8|). In fact, we shall see below that 
for a linear power-spectrum with n < the saddle-point 6l 
is not the global minimum of the action. Actually, in this 
case the action is no longer bounded from below. Then, 
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the steepest-descent method described above is a priori no 
longer justified. In fact, a specific study shows that it is 
still useful but it requires some care. We shall come back 
to this point in the next section. 



the other hand, for j/s < y < we see that we have two 
intersection points t_ < r, < r-|_ < 0. Thus the function 
T(y) is bivaluate over this range. 



3.5. Geometrical construction 

In order to get an intuitive picture of the generating func- 
tion defined by the system (^^ it is convenient to devise 
a geometrical construction which yields T{y) and ip{y). 
First, we note that the first line of eq.JTOl) simply states 
that the implicit function t(i/) is given by the intersection 
of the straight line r/y of variable slope l/y with the fixed 
generating function —G'{t). 




Fig. 2. Construction of the function r(y) for a linear 
power-spectrum with n = — 1 in a critical density uni- 
verse. To each value of y is associated the abscissa r of 
the intersection of the straight line r/y with the curve 
—G'{t). The dashed-line with the triangle shows the loca- 
tion of the singularity Ts — —1.4, ys — —0.44. 

This construction is shown in Fig.|| for the case of a 
linear power-spectrum with n = — 1 in a critical density 
universe. From eq.(|7l|) we obtain the inverse t{Q) as: 



Tig) = -ii + g) 



(n+3)/6 



{i + T)-'[i + g] 



(75) 



where (1+JF)^^ is the inverse of the function 1 + JF. Then, 
using eq.(^) and eq.(^), keeping in mind that g is the 
actual density contrast (eq.JG^)), we get from eq.(|75|) the 
asymptotic behaviours: 

high densities : r ^ -oo, ^ ^ oo : ^ - (_^)6/(«+3) (-^g^ 



and: 



low densities : r — + cx3, ^ ^ — 1 : (1+^) ~ t 



-6/(l-n) 



(77) 



where we did not write positive numerical multiplicative 
factors of order unity. Thus, for large negative r we have 
—g' ^ (— r)'^^~")/'^^+"^ which increases faster than \t\ for 
n < 0. This implies that there is a minimum value j/s < 
of y for which the straight line r/y can intersect the curve 
— ^/'(t). This is shown by the dashed-line in Fig.||. On 
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Fig. 3. The function r(j/) for a linear power-spectrum with 
n = — 1 in a critical density universe. The triangle shows 
the location of the singularity Tg ~ —1.4, ys — —0.44. 

We display in Fig.ljthe function T{y) we obtain for the 
case shown in Fig.||. Note that the function 2/(t) is well- 
behaved and shows no singularity. The singularity (r^ , j/s ) 
is given by the point where dr/dy = oo. From eq.([70[) this 
condition also reads: 



S'(t.) = T,g"(T,), y, = 



1 



g'iTs) g"iTs)' 



(78) 



The advantage of the geometrical construction dis- 
played in Fig.g is that it shows at once the location of 
possible singular points. However, there is an alternative 
geometrical construction which directly yields the value 
of the generating function ip{y). Indeed, we can also write 
the full system ^u) as the one equation: 



(p{y) — TmTi[Sy{T)] with Sy{T) 



yGir) + Y- 



(79) 



Of course, the minimum which appears in eq.(|79|) ex- 
presses the fact that the saddle-point we look for is the 
minimum of the action S'[(5l]. Note that eq.([79|) means 
that we minimize the action S over the subspace of spher- 
ical linear states S^ of the form (pTl) parameterized by r 
(or Sl,Rl). However, we must actually minimize the ac- 
tion over all possible states 5^. The rigorous justification 
of eq.(79) was obtained in Sect. 3. 2 and Sect. 3.4 where we 
showed that the spherical minimum described by eq.([79|) 
is also a minimum with respect to transverse directions. 
Then, we can give the following geometrical solution of 
eq-dzi). 

First, for y > we note that the point r where Sy{T) is 
minimum is also the point where Sy{T)/y = g{T)+T^/{2y) 
is minimum. Then, this point is simply given by the first 
contact of the parabola h — r^/(2y), of varying height h, 
with the curve 5(t), starting from below at h — — oo. 
Then, the minimum of the action is given by ip{y) = y x h 
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Fig. 4. Construction of the function ^(y) for a linear 
power-spectrum with n = — 1 in a critical density uni- 
verse. The quantity ip{y)ly is simply the height h of the 
parabola h — T^/(2y) at the first contact with the curve 
Q{t). For y > (i.e. underdensities) we start from below 
at /i = — CK) while for y < (i.e. overdensities) we start 
from above at /i = -l-oo. 



at this point. Second, for y < we need the maximum of 
Sy{T)/y. This is given by the first contact of the parabola 
h — T^/(2j/) with the curve Q{t), starting from above at 
h — +00, and we have again Lp{y) = y x h. This construc- 
tion is displayed in Fig.H. In particular, it is clear that for 
small y the parabola are very narrow and we get only one 
contact point at r ~ y as we probe the small-r part of the 
curve G{t) where G{t) ~ — r. That is the curvature of the 
parabola gets very large with y — > while the curvature of 
G{t) is finite. This is the essence of the discussion below 



eq.(73). Thus, this geometrical construction gives at once 
the value of the generating function ip{y). In particular, 
one can see at a glance from the curve G{t) the behaviour 
of (p{y). Note that if there exists a singular point (r^, y^), 
as in Fig.y, the minimum obtained for small negative y is 
only a local minimum. We shall come back to this point 
in Sect.|j[ 

Finally, we display in Fig.|| the generating function 
(p(y) obtained for this same case n = — 1 (and fim = 
1). The feature at y = ys shows that (p{y) is singular at 
this point. The curve drawn in Fig|5|is obtained from the 
parametric system (70) with — oo < t < c». As explained 
from Fig.||the function T{y) is bivaluate over y^ < y < 0. 
This also applies to ip{y). The branch which runs through 
the origin in Fig.g| corresponds to r > r^ while the upper 
branch over ys < y < which starts almost vertically at 
ys corresponds to r < r^.. 



3.6. Calculation of the pdfV{6R) 



The generating function (p{y) was obtained in Sect.|3.4|for 
real values of y, using a steepest-descent method. Then, 



Fig. 5. The generating function ip{y) for a linear power- 
spectrum with ri = —1 in a critical density universe. The 
feature at y = y^ — —0.44 clearly displays the singularity 
of (p{y) at this point. This "regularized" generating func- 
tion ip{y) exhibits a regular branch at the origin and it is 
bivaluate over the range ys < y < 0. 



using eq.(|2l|) and eq. (P6|) we obtain the pdiV{5ii) through 
the inverse Laplace transform: 



ViSR) = 



-\-ioo 



dy 



2TTia^{R) 



AvSR->fiiv)]/'y^iR) 



(80) 



Next, in order to compute numerically the pdf V{5fi) 
through eq.(80) we need to continue the function (p(y) 
over the complex plane. This is simply done by using the 
same implicit system (70) for complex values of y. Note 
that the function G{t) is analytic over the region of in- 
terest. Indeed, the function T{Sl) defined in eq.(p9|) and 
eq. (^) is analytic, with singular points such that cos 9 = 1 
or cosh?7 = 1 (with 6^0 and rj ^Q). However, the singu- 
larity at 6l — 1-68 (i.e. for 9 = 27r) is repelled to t = —oo, 
as shown in eq.([75|) and eq.([76|). Hence the function G{t) 
has no singularities along the real axis. 

Then, we need to specify the integration path over y 
in eq.(pO[). It intersects the real axis at the saddle-point 
{tcVc) given by: 



(yc) = with xiv) = V 5r~ viv)- 



dy 



From eq.(|7d) we have: 



G{Tc 



5r since ip'{y) = Gir). 



(81) 



(82) 



Thus, we see from eq.(p8) that this triplet {6R,Tc,yc) 
is also the triplet {5ji,T,y) we obtained in Sect.|3.2| and 
Sect. 3.4 to get ip{y). Of course, this is required by self- 
consistency. It simply means that V{Sr) at the point Sr is 
governed by the neighbourhood of the saddle-point Sl (x) 
obtained in Sect. 3^, which obeys T{5l.Rl) = ^R- This 
result agrees with intuitive expectations. Then, the inte- 
gration path in the complex plane is set by the constraint 
Im(x) = with the requirement that Re(x) decreases on 
both sides from its value at yc. This corresponds to the 
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steepest-descent path. The second derivative of the factor 
Xis: 

x"{y) = -v"{y). (83) 

As seen in Fig.ra we have f" {y) < (if there is no singular- 
ity j/s)- In fact, for \y\ ^ 1 we get: (p"{y) ~ —1. Therefore, 
the path of steepest-descent is orthogonal to the real axis 
at the point yc. Moreover, it is symmetric about the real 
axis, which clearly shows from eq.(|80|) that the result for 
V{Sji) is real, since ip{y*) — (p{y)* . Thus, we can write 
eq.(|80|) as: 



V{5b) = Im 



Ay 



^:a^R) 



,[ySn-<p(v)]/a\R) 



(84) 



where we only integrate over the upper half-plane Im(j/) > 
0. Note that in the quasi-linear limit cr — > the contribu- 
tion to the integral (80|) only comes from an infinitesimal 
neighbourhood of the saddle-point (tc, yc) around the real 
axis. Therefore, we could try a steepest-descent approxi- 
mation for eq.(pO|). However, the agreement with the re- 
sults of numerical simulations is better if wc numerically 
compute the exact integra l (pOJ ). Thus, in the following wc 
compute the full integral (|80|). 

Note that this means that it is better to approximate 
the generating function f{y), and next to use the exact 
inverse Laplace transform (|80|), rather than to directly 
approximate the pdf 'P{6ii). This can be understood as 
follows. In the limit cr — > a steepest-descent approxima- 
tion to eq.(pfl) would be fully justified (it is actually exact 
in this limit). However, it is clear that if we use the re- 
sults obtained for P{6ii) in the limit a ^ (that is, we 
assume that we have obtained in some way the behaviour 
of V{5ii) at all points Sr at leading-order in this limit) 
for a finite value cr > we can generically expect that the 
moments (^|j) obtained from this approximate pdf exhibit 
sub-leading terms which are not correct. In particular, this 
means that we would get: (1) = 1 -f o(l) and {6b) = o(l) 
where o(l) stands for a term which vanishes in the limit 
cr — > 0. Generically, we may expect this term to be of order 
unity when a '^ \. This implies that for small but finite a 
the pdf is not exactly normalized to unity and the mean 
{5ii) is not exactly zero. By contrast, using the exact in- 
verse Laplace transform ( |80| ) with the generating function 
^p{y) (obtained in the limit a —t 0) ensures that we have 
for any finite a the exact integrals (1) ~ 1 and {Sr) = 0. 
Thus, the normalization and the mean are always correct. 
This result can be obtained from the expansion (p3) and 



eq.(26) which shows that in order to have the exact mo- 
ments of order and 1 we only need ip{y) to be quadratic 
in y for y — > 0. Of course, this is the case since from eq.([70|) 
we have the expansion tf{y) — —2/^/2 + .... On the other 
hand, eq.(]80|) imphes that {6'j^i = a'^ for any a. 

The procedure we described above allows us to com- 
pute the pdf V{5b) in the quasi- linear regime, using the 
steepest-descent method developed in the previous sec- 
tions. However, when the function Q{t) grows faster than 
r^ as r — > — oo a singularity ys shows up in the generat- 
ing function ^p{y) and matters are slightly more involved. 



First, we note that for such functions Q{t), which corre- 
sponds to n < as shown by eq.(76), there is no global 
minimunr of the action S[5l\ for negative y. This is clear 
from the construction of Fig.0. Indeed, it is obvious that 
the contact point shown in FigBat r ~ —0.5 for the upper 
parabola is only a local minimum and there is no global 
minimum: whatever large h is taken to be, the parabola 
always intersects the curve Q{t). This is also clear from 
eq.d79). Indeed, we now get Sy{T) -^ — oo for r — > — oo. 
Hence the "action" Sy{T) is not bounded from below if 
y < 0. This actually means that the path- integrals (24) 
and (29) diverge for y < 0. Hence the generating functions 
4'{y) and ^{y) exhibit a branch cut along the negative real 
axis. Then, the steepest-descent method described in the 
previous sections must be modified (in fact, there may 
still exist a global minimum if we take into account shell- 
crossing, which appears for large 5l or large negative r, 
but this is irrelevant here). Note that a negative Sr cor- 
responds to positive Tc and y^ as shown by eq.(|82|) and 
Fig.H. Hence this problems only appears when one looks 
for the value of the pdf V{6r) for positive Sr. 



We can note that from a physical point of view the 
pdf V{Sr) at the point 5r should still be governed by 
the saddle-point {tc, Vc) obtained in the previous sections. 
Indeed, it is clear that a non-linear density contrast Sr 
arises from initial conditions close to the spherical saddle- 
point derived in Sect. 3_^ . In fact, there is a straightforward 
trick to show this in a more explicit fashion. Indeed, as we 
explained above the problem is due to the rapid growth 
of the functional Sr[Sl] for large positive Sl (we do not 
consider shell-crossing here since it is not related to this 
problem). Then, instead of looking for the pdf P{Sr) we 
can as well investigate the pdf 'P{5^'') where g is a large 
odd integer. O bvio usly, the steepest-descent method de- 
veloped in Sect. 3.2 can be applied to this new pdf. This 
involves new generating functions ipq{y) = (e~^ « ) and 
fq{y). We again obtain a spherical saddle-point of the 
form (fill) and the implicit system ( |7(]| ) where the new 



function Qq{T) is simply: Qq{T) = Q{Ty^''. Note that the 
saddle-point Sl,Rl associated with a given non-linear den- 
sity contrast Sr does not depend on q. Of course, this was 
to be expected since to a given Sr corresponds a well- 
defined set of initial states (5l(x), whatever we consider 
Sr itself or (5^ ! Then, we see that if we choose a large 
enough value for q the function Gqir) grows more slowly 



than r for r 



-(X). Therefore, we can now apply the 



steepest-descent method as described in the previous sec- 
tion. Note that the new generating function ipq{y) shows 
no singularity ys,g so that we span the whole curve Gqir) 

(hence G{t)). Finally, from ViSjl'^) we can derive V{Sr) 
through a simple change of variables. For all q we obtain 
in this way the same exponential-like cutoff (i.e. the expo- 
nential of a given power of the density) at large densities 
but the multiplicative factor obtained in the limit a —^ 
will usually differ. In other words, in order to get a unique 
and well-defined result we must take into account the de- 
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Fig. 6. The pdf V{Sr) for n = -1, f7„i = 1 and a = 0.74. 
The data p oints (obtained from n umerical simulations) are 
taken from Bernardcau (1994a)[ T he solid line shows the 
theoretical prediction from cq.(|8C|) and eq.([70|). The dot- 
ted line shows the contribution from the unstable saddle- 
point T_ which gives the high-density tail of the pdf, from 
eq.(86). The dashed- line displays the Gaussian with the 
same variance. 



In fact, as described in App.^ and AppJ^ we can di- 
rectly work with the density contrast Sr, even though 
the path-integral ( |24| ) diverges for negative real y. One 
must simply be careful to use appropriate integration con- 
tours in the complex plane when using integrals like (pOJ). 
Thus, the integration paths we use in eq.(pO|) are shown in 
Fig.|A.2|. They depend on the density contrast Sr and they 



run through the spherical saddle-point derived in Sect. 3.2 



In particular, in agreement with the simple procedure de- 
scribed above based on Sj/', for large Sr the pdf is gov- 
erned by the saddle-point {tcHc) given by eq. ([8^) . This 
means that for high density contrasts we span the upper 
branch of (p{y) shown in FigJ5[ Since the action S[S]^] is 
not bounded from below for y^ < y < the two saddle- 
points r_ and r_|_ (with t- < Ts < t^ < 0) obtained in 



Sect. 3.5 are not the global minimum of the action (which 
does not exist). The saddle-point t_|_ (which corresponds 
to the branch of f^y) which runs through the origin in 
Fig.@) is only a local minimum. On the other hand, the 
point T_ (the upper branch of ip{y)) is an unstable saddle- 
point: it is a local maximum of the action. However, as 
shown in AppH and App.H the Laplace transform ^{y) 
and the pdiViSn) are still governed by these saddle-points 
in the quasi-linear regime. In particular, the saddle-point 
T_ yields the high-density tail of the pdf. Note more- 
over that the generating function ip{y) obtained in eq.([70|) 
and shown in Fig.^ is not the exact generating function. 
Indeed, as we noticed above the actual generating function 
shows a branch cut on the real negative axis (i.e. for y < 
and not only y < y^). This is also explained in App.|A|. 



Finally, we show in Fig.O the pdf V{Sii) obtained for 
cr = 0.74 and n = —1 in the quasi- linear limit. The 
solid line is the result obtained from eq.(pn) and eq.([70|). 
That is, we use the branch which runs through the ori- 
gin of the generating function defined by eq.([70[) and 
displayed in FigJ5|. This curve was already obtained in 
Bernardeau (1994a) through a perturbative method (see 
also Sect. 4.1). Since this function exhibits a branch cut 
on the real axis at y < ys it yields an exponential high- 
density tail of the form 1^(6^) ~ e''=''«/'^'(^) (note ys < 0) 
as can be seen from eq.(p^) (the integral is dominated by 
y ~ ys), see also App.^This describes the pdf V{Sr) 
for Sr < Q{ts) in the quasi-linear limit. However, for 
larger density contrasts the pdf is governed by the unsta- 
ble saddle-point r_ , which also led to the upper branch of 
the "regularized" f{y) in Fig.^. Then, as shown in App.^, 
by closing the integration contour onto the negative real 
axis the inverse Laplace transform (pG) can be written as: 



dy 



T^i^n) = -^ ^"""'^ Im[V(2/)] 



(85) 



Here we noted 2i lm['0(y)] the discontinuity of the gener- 
ating function ^{y) = 'ip{y/a^) — e"'^^^^/'^ along the neg- 



ative real axis, see also eq.(A.14) and eq.(A.l(:). We com- 
pute the pdf V{Sr) in App.^ through a steepest-descent 
approximation, which yields: 



V{5r) 



-T^/{2a\R)) 



27ra 1 + 5r \g'{T) 



(86) 



The pdf obtained in this way is shown by the dotted-line 
in Fig.g. It yields a smoother cutoff for the high-density 
tail of the pdf. Indeed, using the asymptotic form ( |76| ) we 
get from eq.(|8C 



(5fl,»l: V{5R)^e- 



s(" + 3)/3 



/a\R) 



(87) 



where we did not write positive multiplicative factors of 
order unity in the exponent. Thus, we see that for n < 0, 
where the singular point {ts, ys) appears, the pdf exhibits 
a high-density cutoff which is shallower than a pure ex- 
ponential. Of course, this clearly implies that the Laplace 
transform Tp{y) defined in eq.(po|) shows a branch cut on 
the real negative axis since the integral in the r.h.s. of 
cq.(|20|) must diverge for ?/ < 0. This is at the origin of 
the additional difficulties encountered in the case n < 0, 
where an artificial singularity ys appears. This is discussed 
in App.^. The form ( |87| ) is actually valid for all n (i.e. 
also for n > where the analysis is simpler since the rel- 
evant saddle-point t is really the global minimum of the 
action) . It explicitly shows that the high-density cutoff of 
the pdf depends on the initial conditions (the slope n of 
the linear power-spectrum). This is actually quite natu- 
r al and it agrees with a simple intuitive spherical model 
( Valageas (1998) ), as discussed below in Sect.^. We note 
in Fig.y that for a = 0.74 the unstable saddle-point con- 
tribution only dominates V{Sr) for Sr^ 8 which is larger 
than g{Ts) — 3. This is due to the finite value of cr. Thus, in 
the limit cr — > 0, as soon as Sr > QIts) this contribution is 
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larger than the result one would get by using for ^p{y) only 
the branch which runs through the origin in Fig.H. Here, 
we must note that for large density contrasts the form 
( pq ) is no longer valid since shell-crossing comes into play. 
Then, one must take into account virialization processes. 
However, we shall not study these very high densities here 
since we restrict ourselves to the quasi-linear regime where 
such events are extremely rare. Finally, we can see in Fig.ra 
that the quasi-linear limit provides a very good estimate 
of V{5r) up to ct '--^ 1, which is a rather large value. Note 
also the strong departure of the pdf from the Gaussian 
(shown as the dashed-line). 



4. Comparison with previous results 

4.1. Perturbative methods 

Eventually, we point out that the results we obtained in 
Sect.^ partly agree with the standard results derived from 
a perturbative hydrodynamical approach. Indeed, the sys- 
tem dTG) which gives the generating function 93(2/) was 
also obtained by Bernardeau (1994a). This result was de- 
rived from a perturbative expansion of the density field 
over the linear growing mode, substituted into the equa- 
tions of motion of the hydrodynamical description. Note 
that our calculation does not involve the hydrodynamical 
approximation: it is based on the coUisionless Boltzmann 



equation. However, as explained in paper ] the perturba- 



tive expansions obtained in both approaches actually co- 
incide. Hence it is not surprising that we recover most of 
the results of Bernardeau (1994a)[ 

On the other hand, we stress that the method we pre- 
sented in this article is actually non-perturbative. In par- 
ticular, we did not need to assume that the density field 
can be written as a pertur bative ex pans ion. This is impor- 
tant since as explained in paper I and paper V this per- 
turbative expansion actually diverges (it is only asymp- 
totic). Moreover, our calculation directly provides the pdf 
V{Sfi) in Eulerian space and we do not need to go from 
Lagrangian space to Eulerian space, which is a delicate 
step in the usual method. In particular, there is no need to 
apply any "smoothing" a posteriori: we directly obtain the 
pdf of the density field at a given scale R which enters into 
the formulation of the problem itself. Thus, our calculation 
provides the needed justification of these previous results. 
For instance, if the spherical saddle-point we obtained in 
Sect. 3.2 were only a local minimum of the action and there 
were another deeper minimum for a non-spherical density 
field S']^ of the same order this would show up in our for- 
mulation and we could take into account the contribution 
of this second minimum. By contrast, the perturbative ap- 
proach would not provide this second minimum (the only 
hint of its existence would be that the perturbative series 
diverges, but this is the case anyway for other reasons). 
Fortunately, as shown in Sect. 3^ matters are simpler than 
this and in the quasi-linear regime the pdiV{6H) is indeed 
governed by this trivial spherical saddle-point. 



Note however that these previous works based on the 
perturbative approach always used the implicit system 
( [70| ) for the generating function. More precisely, they used 
the branch of the generating function ip{y) which runs 
through the origin in Fig||. As discussed in Sect. 3.6 this 
means that for n < they get a mere exponential cutoff 
V{Sr) ~ e^=''^/°'^(-'^^. As shown in App.^ and discussed in 
Sect. 3.6 this is actually incorrect. Indeed, the high-density 
cutoff of the pdf is of the form (p^) (until shell-crossing 
occurs) and the actual generating function (p{y) shows a 
branch cut along the whole real negative axis (and not 
only for y < ys) if we disregard shell-crossing. Thus, we 
see that for n < the perturbative approach fails beyond 
Sr > G(ts) and the high-density tail actually requires a 
non-perturbative treatment. In fact, the "resummation" 
of the perturbative theory at leading order performed in 
Bernardeau (1992) yields the implicit system ( |7C| ) which 
remembers the existence of the non-perturbative unsta- 
ble saddle-point t_ < Ts- Indeed, eq.([70|) can also be ex- 
tended to r < Ts where it yields the upper branch of (p{y). 
However, in order to use the information contained in this 
upper branch one needs the non-perturbative method de- 
scribed in this article, which provides the integration con- 
tour required to take into account the contribution of this 
unstable saddle-point, see Fig, A. 2, and which gives the 
full justification of this procedure. 

This unstable saddle-point modifies the pdiViSu) for 
^R > G{'''s)- This implies that the moments and the cu- 
mulants of the pdf are also changed. Therefore, at finite 
(T they are not given by the expansion around y = as in 
eq. (|28|) of the "regularized" generatin g fu nction (p{y) ob- 
tained from ( |70|) (as discussed in Sect.3_^ the exact gener- 
ating function if{y) for finite a is not regular at the origin). 
Moreover, the additional contribution to the moments of 
the pdf which arises from this shallower high-density cut- 
off is non-perturbative. Indeed, from eq.(87) we see that 
the change induced by this correction is of order: 



M^'r) 



,-e(r,)(" + 3)/V<T^(ii) 



(88) 



since the pdf is only modified for Sr > G{ts)- It is clear 
that this correction cannot be obtained by a direct per- 
turbative treatment, where one would derive the moments 
of the pdf by computing in serial order the terms which 
arise from a perturbative expansion of the density field. 

Indeed, the correction 



as m 



Goroff ct al. (1986) 



j) can- 
not be written as a perturbative expansion over pow- 
ers of a (it vanishes faster than any power of a in the 
limit a —^ 0). This is why the results of this iterative 
perturbative method actually coincide with the the re- 
sults obtained from eq.(Eq) by expanding the "regular- 
ized" generating function (p{y) defined by eq.lTQ). Thus, 
we see that a standard iterative p erturbative method, 

cannot give the pdf 
In particular, these 



(1986) 



as described in [Goroff et al. 

V{dR) for Sr > g{T,) with n < 0. 

features imply that the perturbative series diverge. Note 

that these arguments do not involve shell-crossing. Thus, 

for n < the divergence of the perturbative series is not 
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due to shell-crossing which comes into play at larger den- 
sities. Of course, as noticed in paper I and paper V, for 
n > shell-crossing also leads to a divergence of the per- 
turbative series. Then, we must point out that for large 
density contrasts {Sn ^ 100) shell-crossing must be taken 
into account as the spherical collapse solution is no longer 
described by eq.(B3). Thus, the results obtained in this 
article only apply to smaller density contrasts. 

4.2. Spherical model 

We can also note that the saddle-point method developed 
in this article is somewhat similar to the spherical model 



presented in Valageas (1998) (§2). This model is based on 
the "educated guess" : 



dS{l + S)V{S) 



dSLVLiSL). 



(89) 



Sl.r^ 



It merely states that the fraction of matter enclosed within 
spherical cells of radius R and density contrast larger 
than Sji (in the actual non-linear density field) is approx- 
imately equal to the fraction of matter which was orig- 
inally enclosed within spherical cells of radius Rl and 
overall linear density contrast larger than 5l,Rl- Here 
Rl and 5l,Ri^ are related to R and 5r as in eq.(p^, 
using the spherical collapse solution. Note that this is 
very close in spirit to the Press-Schechter prescription 
used to estimate the mass func tion of just-collapsed ob- 
ject s ([Press fc Schechter (1974)[) . Then, from eq.(||) and 
eq.(38) one obtains the expression (|8q ) for the pdf 7^(^fl). 
In terms of the more usual variable u = 5L.Ri^/(y{RL) this 
can also be written as eq.(B.9). Hence, we see that the 
modified exponential tail (^7j) can be understood in very 
simple terms. It is directly related to the Gaussian cutoff 
of the linear density fluctuations and to the slope of their 
power-spectrum. 

In Valageas (1998)| we showed that the generating 
function ^p{y) of the quasi-linear regime, defined by 
eq.(70), could be recovered from eq.(p9|) and the calcula- 



tion involved a saddle-point as in the present calculation. 
However, that previous work was a simple phenomenolog- 
ical study, based on a simplified description of the density 
field. By contrast, the present work is a rigorous study 
based on the exact equations of motion and we deal with 
the exact 3-dimensional density field. In particular, we do 
not require the density field to be spherically symmet- 
ric. We recover the results of the simple spherical model 
because the saddle-point is spherically symmetric and at 
leading order the generating function ^{y) is given by the 
value of the action at this point. However, our results 
should differ when we consider higher-order terms. 

Here we must note that, as described in Sect.^|^ and 
App^, for n < we need the prefactor [Det(Ai,.Mj^)]"^/2 
for the high-density tail of the pdf (i.e. Sr > G{ts))- We 
did not derive this determinant in a rigorous manner hence 
the multiplicative factor which appears in eq. ( p6[) may not 
be exact. In fact we can expect a non-zero correction to 



the approximation we used in App.O because the exact 
problem we investigate here shows some important differ- 
ences with the simplified spherical model (p9). Indeed, this 
latter model only involves ordinary integrals and a one- 
dimensional variable S^^ji,^. By contrast, the formation of 
large-scale structures in the universe involves the infinite- 
dimensional variable 5l(x) which leads to a path- integral 
formalism. Then, we can expect the integrations over the 
fiuctuations around the saddle-point to show some dif- 
ferences between both cases. However, as seen in FigJfl 
the expression (pq) should provide a reasonably good ap- 
proximation to the exact high-density tail. In particular, 
it should be sufficient for practical purposes. In fact, it 
is probably even sufficient to use the pdf obtained from 
the "regularized" generating function ip{y) (i.e. the lower 
branch in Fig.|^) through eq.(|80|). Note in any case that 
the exponential term obtained in eq. (p6|) is exact, since 
it only depends on the value of the action S[Si,] at the 
spherical saddle-point derived in Sect. 3.2 and not on the 
second-derivative of the action. 

4.3. Velocity field 

We can note that using perturbative methods as in 



Bernardeau (1994a) or th e approximate spherical model 
(^9[), see Valageas (1998) , it is also possible to derive the 
pdf (and the associated generating function) of the mean 
divergence 9 = (V.v)/d of the peculiar velocity field v 
within spherical cells. We shall not compute explicitly this 
pdf 1^(9) here, using the saddle-point method we devel- 
oped in the previous sections. Indeed, it is clear that we 
must recover the results of the hydrodynamical perturba- 
tive approach (i.e. the same generating function Lpg{y)). 
In fact, as long as the test-field :/(x, p, t) which enters the 
functional Z[i] defined in eq.(^ is spherically symmet- 
ric we can look for a spherical saddle-point. Then, since 
the physics involved is the same as the one which governs 
the behaviour of the pdf V{Sji) we shall recover the same 
spherical saddle-point and the results of the hydrodynam- 
ical perturbative method, with the appropriate modifica- 
tion of the tail arising from large densities, as in eq.(p7|). 
Note that for the divergence 6 the pdf shows an expo- 
nential tail for n = —1 (e.g., Valageas (1998)) so that the 
feature which appeared for n < (i.e., the singularity ys) 
for V{5r) is now obtained for n < — 1 for V{9). 

5. Conclusion 

In this article, we have developed a non-perturbative 
method to obtain the pdf 'P{5fs) of the density contrast 
within spherical cells in the quasi-linear regime. This cor- 
responds to a rare-event limit: the rms fiuctuation a van- 
ishes while the density contrast is kept fixed. Then, a 
saddle-point approximation yields asymptotically exact 
results in this limit. Note that our approach does not 
rely on the hydrodynamical approximation for the equa- 
tions of motion. It is fully consistent with the collisionless 
Boltzmann equation. However, it happens that the spher- 
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ical saddle-point which governs the quasi-hnear regime is 
an exact solution of both formalisms (hydrodynamics and 
Boltzmann equation) . This makes the problem rather sim- 
ple and it does not introduce any approximation. This is 
also consistent with the fact that the perturbative series 
obtained from the hydrodynamical and the Boltzmann 



frameworks are identical, see paper I. Although the nu- 
merical examples described in this article were obtained 
for a critical density universe our method applies to any 
cosmological model. One simply needs to use the relevant 
spherical collapse solution .F(Jl) associated with the re- 
quired values of the cosmological parameters f2ni and Q\. 

Thus, we have recovered most of the results obtained 
by the usual perturbative method for Gaussian primordial 
density fluctuations. This provides a rigorous justification 
of these results. Moreover, we have corrected an error in- 
troduced in these previous works for the high-density tail 
of the pdf for power-spectra with n < 0. This clearly shows 
that one should not ask too much from perturbative meth- 
ods, especially since all perturbative series actually diverge 
which gives room for strong non-perturbative corrections. 

Note that our approach is actually much more intuitive 
and simpler than the perturbative method. In particular, 
the spherical collapse solution of the dynamics appears 
naturally in this framework as a saddle-point of the action, 
simply through the spherical symmetry of the problem. 
This symmetry is due to the homogeneity and isotropy of 
the primordial density fluctuations and to the fact that 
we consider the density contrast Sji within spherical cells. 
Then, we have described a geometrical construction of the 
generating function (p{y) (related to the Laplace transform 
of the pdf) which allows one to see at a glance its main 
features. 

To conclude, we note that the approach developed in 
this article presents the advantage to introduce a method 
which is of standard use in physics. In particular, it makes 
the physics involved rather transparent. Finally, another 
advantage of our approach is that in principle it can 
also be applied to non-Gaussian primordial density fluc- 
tuations. This will be described in a companion paper 
( paper III ). Besides, since it is non-perturbative and it 
does not rely on the hydrodynamical description it could 
also be applied to the non-linear regime. In this case, it 
would give the tails of the pdf 'P{Sr) (the saddle-point 
approximation only yields asymptotic results). We shall 
present a study of this non-linear regime in a future work, 
see paper IV. 



Appendix A: A worked example: the lognormal 
pdf 

Here we apply the steepest-descent method to the log- 
normal probability distribution function. This allows us 
to illustrate on a simple example the features implied by 
pdfs with a rare-event tail which decreases more slowly 
than an exponential cutoff. This also corresponds to gen- 
erating functions ip^y) and (pijj) which exhibit a branch 
cut on the negative real axis. 



In order to facilitate the comparison with the problems 
dealt with in Sect .0 we shall use the same notations as far 
as possible. Thus, from a Gaussian variable 6l with the 
pdf: 



r{SL) 



1 



-5i/(2a=) 



27rCT 



we define the new variable p: 
1 1 



p = e '^ hence V{p) = 



,-(lnp)V(2^^) 



P V27rcr 



(A.l) 



(A.2) 



Thus, Sl is the analog of the linear density contrast while 
p — 1 + S is the non-linear overdensity. In particular, Sl 
runs from — oo to -|-oo while p > 0. Here it will be more 
convenient to study p rather than 5, so as to avoid un- 
necessary factors —1. The main difference with the prob- 
lems studied in SectJS is that S^ and p are simple random 
variables and no longer random fields. Hence the path- 
integrals of the main text are replaced here by ordinary 
integrals. This simplifies the discussion and it will allow us 
to compare the predictions of the steepest-descent method 



with the exa ct p df (A.2). 

The pdf (A^) is a lognormal law. We did not shift the 
mean of the Gaussian variable Sl in order to ensure that 
(p) = 1 since it is irrelevant for our illustrative purposes. 
Moreover, we do have (p) = 1 at the leading order in the 
limit cr — > 0. The mome nts of the pdf V{p) can be easily 
computed from eq.(A.l) which yields: 



{p' 



(A.3) 



As in eq.(M) we can define the generating function: 

/•oo 

i;{y) = (e-yp) ^ dp e"^" r{p). (A.4) 

Jo 



Making the change of variable p -^ Sl we get: 
1 



V'(y) = 



dSr 



,-ye'i'~5i/{2a-') 



2Tra 



(A.5) 



This expression clearly shows that ip{y) is not well-defined 



for negative real y since in this case the integral in eq.(A.5) 



diverges for large positive Sl- This also implies that the 
series expansion of tpiv) diverges for all non-zero y: 



V'(y) 



E 

g=0 



^(.^) 



j^tvl 



q=0 



q 



(A.6) 



Next, let us define as in eq.(Eq) the rescaled generating 
functions (p{y) and ip{y)- 



V'(y) = '0(ycr^ 



'Viyo- )/o- 



(A.7) 



As in eq.(|6q) and eq.(|8|) we also introduce the variable r 
and the function G{t) by: 



= -dL, g{T) = 



e-^ = p. 



Then, we can write eq.(A.5) as: 
dr 



V'(y) 



^[yg(r)+T'/2]/a' 



27rcr 



(A.; 



(A.9) 
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This expression is the analog of the path- integral (E9|) . In 
particular, a steepest-descent method yields again exact 
results in the limit cr ^ 0. Moreover, the saddle-point and 
the value of the generating function ip{y) are again given 
by eq.(|70|). Then, from (p{y) one obtains the pdf V{p) (in 
the quasi-linear limit) through the inverse Laplace trans- 
form Q. 

This steepest-descent method is fully justified for pos- 
itive real y where the integral ( |A.9D converges. Since neg- 
ative S (i.e. p < 1) corresponds to positive Tc and j/c 
this method yields the pdf 7'(/o) for p < 1. However, 
for negative real y the integral diverges. Hence one can- 
not directly apply this procedure for p > 1 since the 
saddle-point yc which would appear in the computation 
of eq.(|80|) would be negative. Nevertheless, the steepest- 
descent method is still useful but it must be applied with 
some care, as we shall describe below. A similar prob- 
lem arises in usual Quantum Field Theory when one tries 
to derive non-perturbative results from path- integrals. 
This leads to the s o-called "instanton" contributions, see 
Zinn- Justin (1989) . However, since some features are spe- 
cific to our case (e.g., the saddle-points Tc are not fixed) 
we shall detail the procedure required by the problem we 
investigate. 

First, we need to perform the analytic continuation of 
ipiy) over the complex plane, starting from real positive 
y. To do so, we must deform the integration path C in 
eq.(A.9) as we change the argument of y so that Re(S') 
remains positive for r ~f — oo, where the "action" S is: 



S^yg{T) + ^=ye-^ + ^. 



(A.IO) 



Since this deformation must be continuous as we increase 
(or decrease) Arg(y) from the case Arg(y) = (where C 
is the real axis) one can easily check that the contour C is 
of the form shown in Fig. A.l. It obeys: 



Re(T) -^ — oo : Im(T) = Arg(y). 



(A.ll) 



For Re(r) -^ +oo we can keep Im(T) — (actually, we 
only need |Arg(T)| < 7r/4). Thus, for y = — jyl -I- iO we 
have the contour C+ with Im(T) = tt for large negative 
Re(r) while for y — —\y\ — iO we have the contour C_ with 
Im(T) — —IT. Note that tp{y + iO) ^ ip{y ~ iO), where y is 
real negative. 

We also display in Fig. A.l the points t_, r^ and t+ 



obtained for a small negative y. As in Sectj3.5|, the point 
T_|_ is a local minimum of the action S along the real axis. 
In fact, it is the global minimum along the integration 
path C. On the other hand, r_ is a local maximum along 
the real axis. This implies that the steepest-descent path 
runs through t+ along the real axis while it runs through 
r_ perpendicularly to the real axis (so that r_ is a local 
minimum). From eq.([78|) the singular point (r^, yg) is given 
by: 
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Fig. A.l. The integration path C over r in the complex 
plane. The paths C+ and C_ shown in the figure corre- 
spond to Arg(?/) = TT and Arg(j/) — — tt. They run from 
Re(T) — — oo (where Im(T) ~ ±7r) up to Re(T) = +oo 
(where Im(T) = 0). They give the integration path for the 
generating function tpiv)- The dashed-line shows the con- 
tour d which yields the discontinuity of ■(/'(y) along the 
real negative axis (branch cut). The saddle-points r_ and 
T+ on the real axis correspond to small y: y^ < y < 0. 



Let us first consider for small y < the sum of both 
contributions: 

2 Re[^(y)] = ^(y + lO) + Jj{y - zO) 
dr 



[y g(r)+r' /2\/o 



C++C- 



27rcr 



(A.13) 



where we note "(/"(y) = ^/'(y -I- iO) for y < 0. Then, we 
see from Fig. A.l that the saddle-point t+ contributes to 
Ke[ip{y)] since both contours C+ and C_ run through this 
point along the real axis in the same direction (their con- 
tributions are equal and they sum up). On the contrary, 
the saddle-point t_ gives no contribution. Indeed, even 
if we deform the contours C+ and C_ so that they run 
vertically through r_ in Fig. [A.l| their contributions are of 
opposite sign (both contours are symmetric with respect 
to the real axis). Next, the discontinuity of 4'{y) along 
the branch cut (i.e. the real negative axis) is given by the 
difference between both integrals: 



2i lm[ip{y)] = ip{y + iO) - ijj{y - iO) 



dr 



-[ye(r)+rV2]/<T^ 



27rcr 



(A.14) 



Indeed, the contour C+ — C_ can be deformed into the 



^1, Vs 



-1/e. 



(A.12) 



contour Ci shown by the dashed line in Fig. A.l. Thus, 
the discontinuity is governed by the saddle-point t_ and 
the contribution from t+ cancels. Note however that this 
behaviour only applies to small negative y, y^ < y < 0, 
where the two real saddle-points t_ and t+ exist. For y = 
ys these two saddle-points merge and for y < y^ we have 
two complex conjugate saddle-points. However, we shall 
not need study this regime. 
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The pdf 7'(p) is obtained from the generating functions 
tpiy) or (p{y) through the inverse Laplace transform as in 
eq.(80). This now reads: 

dy 



Vip) 



aPy/°'~ 



2TTia'^ 



i'iy) 



(A.15) 



Since tp{y) has a branch cut along the real negative axis 
the contour Cy bends around the real negative axis and it 
intersects the real axis at y > 0. Th e int egration path Cy 
is shown by the dashed curve in Fig, A. 2 . As we have seen 
above, the saddle-point r+ gives the dominant contribu- 
tion to Re('0) for small y. If we only take into account this 
contribution, we obtain the generating function ip{y) de- 
fined by eq.(|70|). This corresponds to the branch for y > ys 
which runs through the origin in Fig.^. Note that this 
function ip{y) is no longer singular at y = 0. Indeed, as 
shown above the discontinuity of ipiy) along the branch 
cut (which starts at y = 0^) is given by the saddle-point 
r_ and not by t^, see eq.( A.13) and eq.( |A.14 ). This "regu- 
larized" function (p{y) still exhibits a branch cut for y < y^ 
when the saddle-point t+ gives rise to two complex conju- 
gate saddle-points. Of course, since these generating func- 
tions ip{y) and ip{y) are regular at the origin, their expan- 
sion (A. 6) now converges for \y\ < \ys\. Note indeed that 
the moments of the pdf obtained in this way are no longer 
given by eq.( |A.3D since we only keep the leading order of 
the cumulants in the limit ct — > 0. This describes the pdf 
for small positive density contrasts 6 {p ^ 1). Moreover, 



Thus we integrate the discontinuity along the real nega- 
tive axis. Then, in the limit p —^ oo the integral (A. If) is 
dominated by the behaviour of Im[?/'(y)] in the neighbour- 
hood of y ~ 0". As was shown in eq.( A.14), in this regime 
lni['0(y)] is governed by the saddle-point t_. Indeed, for 
y ^ 0~ we have y^ < y < so that the saddle-points r_ 
and r_(- lie on the real axis. Note that the limit p ^ oo is 
not the limit cr —> 0. In other words, in order to obtain the 
high-density tail of the pdf V{p) one must keep a finite 
(even though small) and study the limit of large p. 



, Im(y) 




Re(y) 



one can check by a direct calculation from eq.(A.3) that 
the expansion around y = of the "regularized" function 
(p{y) defined by eq. (|70|) indeed yields the cumulants of the 
pdf 'P{p) at the leading order in cr ^ 0. 

However, the behaviour of the pdf 7-'(p) for large p 
is not governed by r+ but by the saddle-point t_. First, 
let us note that the pdf obtained by this "regularized" 
function (p{y) exhibits an exponential cutoff of the form 



Fig. A. 2. The integration path Cy (dashed curve) over y 
in the complex plane. For p < 1 this integration path can 
be deformed into the contour shown by the solid curve on 
the right which runs through the saddle-point yc on the 
real axis. For 1 < p < p^ or p > ps we use the two contours 
shown by the solid curves on the left (see main text). 

We shall now derive the high-density tail of the pdf 
V{p). First, as explained above we need Im[i/)(y)] for 
y ^ 0^. This is given by the saddle-point r_ in eq.( |A.14 ). 
Thus, a Gaussian integration yields: 



2zIm[^(y)] = -?:[-l-yg"(T_)] 



-1/2 ^-Sy/a' 



(A.17) 



oVsPJa 



. This is obvious from eq.(|A.15D. Indeed, ru \ ■ ■ ■ j , _,„ k „, . , 

, . j„ ,. r ,■: ' I \ J where b^ = oKt-) is given m eq.flA.lO). then, as m the 



since these "regularized" generating functions '^(j)) and 
'(/'(y) show a bran ch cut for y < ys the integration path 
over y in eq. (|A.15 ) is bent around this branch cut and for 
large p (or 5) the integral is dominated by y ~ y^ since 
the other parts of the path with Re(y) < ys become ex- 
ponentially small as e^^'^ with respect to ev^pl" (e.g.. 



Bernardeau (1992)). Second, as we have noticed above the 
exact generating functions ■(/'(y) and (p(y) actually show a 
branch cut along the real negative axis for y < 0. Then, 
it is clear that for large p the pdf is governed by the sin- 
gularity at y = 0. Indeed, formally this leads to a cutoff 
V{p) ~ ey=>Pl° with y^ -^ for large p. This actually cor- 
responds to a pdf with a large density tail which decreases 
more slowly than an exponential. Of course, we can check 

that this agrees with eq.( |A.2| ). 

This property can be derived as follows from eq.( A.15| ). 
If we make the upper and lower branches to get very close 
to the real negative axis we can use eq.(A.14) to write: 



discussion of eq.(piD the pdf V^p) in eq.( |A.16[ ) is given by 
an integration around the saddle-point {yc,Tc). However, 
the integration path now runs along the real negative axis 
and we have ys < yc < and Tc < Tg. It is interesting to 
consider the steepest-descent approximation for this inte- 
gral too. This yields: 



V{p) 



1 



1 



1 



.-rV(2<T^) 



27ra ^-l-yg"{T) ^^"{y) 



(A.18) 



where r obeys Q{t) — p while y and ^p{y) are given by the 
implicit system ([70|). More precisely, (p{y) is now described 
by the upper branch in Fig.0 which runs over ys < y < 0. 
From eq.(|7C|) we also obtain: 



(A.19) 



'p"iy)^G'ir)^. -i-yG"{r) = g'{r)^. 

ay Olt 



This yields: 



V{p) 



_y_ f.pv/'y''- TmrWi, 



ln\[i){y)]. 



(A.16) V{p) 



1 



1 



-rV(2.^) 



2w \g'{T) 



(A.20) 
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Thus, we actually recover the exact pdf (A. 2) since from 
eq.(A.8) we have: r = — Inp and G'{t) — ~p. This is 
natural for large p since the steepest-descent method de- 
scribed above becomes exact in this limit. For smaller p 



the pdf given by eq.(A.15), where we perform the exact 
numerical integration, is no longer equal to the exact pdf 
( [A.2| ). This is the generic case. The pdf obtained by the 
steepest-descent method only yields asymptotically exact 
results in the two limits cr —> or p ^ oo. 



— -2 — 




Fig. A. 3. The pdf 1^(6) for the lognormal case wit h a ~ 
0.74. The dashed-curve shows the exact pdf from eq.(A.2). 
The solid curve shows the results obtained by the steepest- 
descent method. The break at 5 ~ 2 corresponds to the 
transition to the high-density tail which is governed by 
the saddle-point r_. 



Finally, using the steepest-descent method described 
above we can estimate the pdf V{p). As discussed above, 
the exact generating functions ^p{y) and ip{y) show a 
branch cut on th e real negative axis so the integration 
path over y in eq.(A.15) is bent around this cut as shown 
by the dashed curve in Fig. A. 2 , For p < 1 this contour can 
be deformed in the path shown by the right solid curve 
which runs through the saddle-point j/c > on the real 
axis. For p > 1 matters are more intricate as discussed 
above. For 1 < p < p^, where ps = GiTs) = 1/e, we only 
take into account the "regularized" part of (p{y) which is 
described by eq.(^. This function is regular over y > ys 
and we use the integration path shown by the solid curve 
labeled "1 < p < p^" , which runs through a saddle-point 
J/s < 2/c < on the real axis. Next, for large densities 
p > Pa we take into account the exact branch cut along 
the real axis in order to obtain the correct high density 
tail. This is shown by the contour labeled "p > ps". We 
split this path into two parts. The first one around the 
real negative axis with j/s < 2/ < is governed by the 
saddle-point t_ which leads to a saddle-point j/c which 
goes to 0~ for large p. The second part over the range 
— oo < Re(y) < ys is computed using eq.(70) for (f{y). 
Then, we simply take the largest of these two contribu- 
tions to estimate 'P(p). In fact, for p > Ps the first part is 



the largest one in the limit cr — > 0, as expected. For ct ^ 
or p — > oo it is actually infinitely larger. 

We display in Fig. |A.3| our results for the pdf 7'((5) with 
(7 = 0.74 (we show V{5) rather than V{p) in order to com- 
pare with the main section). The results obtained by the 
steepest-descent method are shown by the solid curve. The 
knee at (5 ~ 2 corresponds to the transition from the "reg- 
ularized" (/3(j/) to the exact branch cut at j/ < 0. Thus, for 
(5^2 the pdf is governed by the saddle-point t_ and the 
neighbourhood of y = 0~. We can check in the figure that 
this estimate is indeed exact in the limit p — > oo, as we 
proved in eq.(A.20). For smaller p the pdf is governed by 
the saddle-point t_|_, that is Lp{y) is described by the branch 
which runs through the origin given by eq.(|70|). Thus, we 
see that the steepest-descent method provides reasonably 
good results up to cr ~ 0.74. However, the agreement is 
not as good as in Fig.o for the actual pdf V{Sr) which 
arises from gravitational clustering. 



Appendix B: Application to path-integrals 

We have described in App.^ how to apply the steepest- 
descent method in the case where the function G{t) grows 



faster than t'^ for r 



-oo. As noticed in SectJS.q and 



Sect. 3.6, this corresponds to linear power-spectra P{k) 
with n < for the problem of gravitational clustering 
which we investigate in this article. The problem stud- 
ied in App.H actually involved ordinary one-dimensional 
integrals but the arguments can be generalized to path- 
integrals. Note that in our case the function G{t) does 
not grow as e~'^ (as in App.|A[) but as a power-law, see 
eq. (|7q) . Then, the contour C in the complex plane over r 
which was shown in FigjA.l] is now given by: 



Re(T) 



-oo 



Arg(T) 



-Arg(y). 



(B.l) 



However, this does not change the behaviour described in 
App.0 



Thus, we can directly apply to the path-integral (29) 
the procedure detailed in AppH. However, for n < 
where the high-density tail is governed by the saddle- 
point r_, t he transposition of the Gaussian integration 
of eq.( A.14 ) which yielded eq.( A.17 ) now gives a factor: 



2i Im[i/'(y)] = (DetA^ 



,1/2 



(DetM^y) 



-1/2 -Sy/a^{R) 



(B.2) 



as in eq. (|63|) . Note that since the spherical saddle-point 
(5l(x) (i.e. T_) is now a local maximum of the action S[Sl] 
the determinant DetMj, is negative so that the square root 

1/2 

(DetMj,) ' yields a factor i. Thus, the numerical factor 



which appeared in eq.(A.17) is replaced by a determinant: 



l + yG"{T^)^'D with V = r>ct{AL.My 



(B.3) 



In order to compute the pdf V{5ii) we only need to eval- 
uate this determinant T). Indeed, the minimum Sy of the 
action is still given by eq. (|7C|) where we use the upper 
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branch defined over the range ys < y < shown in Fig.ra. 
From eq.(p3h we have: 



(AL.My)(xi,X2) =(5d(xi -X2) 



CT^R) 



/ dx Al(xi,x) 



SH^r) 



5(<5L(x))5(5i(x2)) 



(B.4) 



Thus, we need the second-order derivative of the func- 
tional (Jfl^f^L] taken at the spherical saddle-point charac- 
terized by 5l,Rl (or ■'')• Since we do not know the explicit 
form of the functional 5]:([5l\ we shall only obtain a simple 
estimate of the determinant T). However, as we can see in 
Fig.g the behaviour of the discontinuity of the generating 
functions i/'Cy) ^-iid ^p{y) near y = Q~ only dominates the 
pdf V{Sji) for rather large density contrasts {S ^ 8) for 
a = 0.74. Hence we do not need the high-density tail of 
the pdf with a high accuracy if we study the quasi-linear 
regime. 

In order to estimate the determinant V we could try to 
use second-order pertur batio n theory. Indeed, the second- 
order derivative in cq.(B^) becomes negligible on very 
large scales x ^ Rl, X2 ^ Rl- Moreover, the radial profile 
of the spherical saddle-point is almost flat in the inner re- 
gion R' < Rl ■ Hence we might estimate this second-order 
derivative, which should be taken at the point Sl{x) given 
by eg. (|6l|) , by its value at the point (5l(x) = 5l,Rl (i-S- 
constant density contrast). Then, we simply need to in- 
vestigate the second-order perturbation theory in a back- 
ground universe characterized by a higher mean density: 
p — > p(l -I- Sr). Here 5r = g{T) = T{6l,Rl) is the ac- 
tual non-linear density contrast of the spherical saddle- 
point. Unfortunately, this procedure cannot give mean- 
ingful results. Indeed, it is well known that perturbation 
theory leads to divergent quantities when one goes beyond 
leading order terms (e.g., Scoccimarro fc Frieman (1996)[ ). 
Then, it is easy to check that the calculation of the de- 
terminant D from eq.( |B.4[ ) with the use of the second- 
order term ^(^^(x) for the density field S{x) (written as 
an expansion over <5l(x)) gives rise to such divergences. 
In fact, as shown in [paper V we can check that using a 
perturbative approach to evaluate the fluctuations of the 
action S[6l] around the saddle-point we exactly recover 
the divergences obtained from standard direct perturba- 
tive methods. 

As a consequence, we shall use the following approxi- 
mation for the generating function ip{y). By analogy with 
the case of ordinary integrals studied in App.H we replace 
in a first step the determinant I? by a factor 1 + yQ"{T^), 
see eq.(B.5). This takes into account the dependence of 
the non-linear density contrast on the local linear density 
contrast Sl,Rl- However, a new physical process which 
did not appear in App.H occurs in the context of cos- 
mology: the expansion of the background universe. This 
leads to a dilution of the high-density tail of the prob- 
ability distribution V{5r). Indeed, let us consider for a 
moment the following local model. At a time ii we mark 
the comoving coordinates x where the linear density con- 
trast 5i,fl^^i(x) over the cell V^i centered on x is above 



some threshold 5ci , which corresponds to a non-linear den- 
sity contrast Ad = !F{5ci)- This fills a fraction Fi of 
the volume of the universe. At a later time ^2, the same 
fraction Fi of the universe now shows density contrasts 
above 5c2 — 5ciD+[t2) / D+{ti) and Ac2 = ^{Sc2), where 
D^{t) is the linear growing mode. However, it happens 
that in fact the regions over the non-linear threshold Ac2 
no longer fill a fraction Fi of the universe. Indeed, while 
their density increases these regions also depart from the 
mean background expansion and they actually contract 
in comoving coordinates (for positive density contrasts). 
Thus, we have F2 = Fi(i?2/i?i)^ = Fi(l + Aei)/(l + Ae2). 
Therefore, we add a dilution factor 1/(1 + 5r) to the gen- 
erating function ip(y)- Hence we write: 



2i Im[i^(y)] 



1 



-Sy/a\R) 



l + '5fl y/T+Wir) 



(B.5) 



where 5r, = Q{t) and T(y) is given by the upper branch in 
Fig.^. Next, the pdf P{6r) defined by t he inv erse Laplace 
transform ( pi| ) can be written as in eq.(A.16) in terms of 
Tpiy). This yields: 



VidR) 



1 



1 



" dy 



.[ySn-Sy]/i7^ 



.(B.6) 



Then, in the limit cr — )■ we can evaluate this integral 
by an ordinary steepest-descent method, as in eq.(A.18). 
Thus we obtain: 



V{Sr) 



1 



1 



1 



-rV(2a^(fl)) 



27ra 1 + 6r |e'(r)| 



(B.7) 



The expression (B.7) can actually be recovered from a very 
simple spherical model detailed in Valageas (1998) . Thus, 
let us define the variable v which describes linear density 
fluctuation by: 



OLMl 

<y{RL 



a{R)- 



Then, the expression (B.7) can also be written: 
1 1 di/ 



ViSR) 



-1.72 



/27r 1 + (5fl d6R 



(B.i 



(B.9) 



where we used: G'{t) = dSR/dr. The relation (B.9) 
agrees exactly with the spherical model studied in 



Valageas (1998), which we discuss here in Sect. 4.2 



This model is similar to the Press-Schechter approach 
( Press &: Schechter (1974 )| ) used to estimate the mass 
function of collapsed halos (without the factor 2). Then, 
eq.(B.9) should be sufficient for practical purposes, as seen 
in FigJ^. Note that although we did not obta in a rigorous 
derivation of the multiplicative factor in cq.(B.£) the ex- 
ponential cutoff is exact. 
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